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Abstract 

We put forward the existence of a rich variety of fully stationary vortex structures, termed 
H-clusters, made of an increasing number of vortices nested in paraxial wave fields confined by 
trapping potentials. However, we show that the constituent vortices are globally linked, rather 
than products of independent vortices. Also, they always feature a monopolar global wave front 
and exist in nonlinear systems, such as Bose-Einstein condensates. Clusters with multipolar global 
wave fronts are non-stationary or at best flipping. 

PACS numbers: 42.25.-p 
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Singular wave structures, that contain topological wave front dislocations are ubiqui- 
tous in many branches of classical and quantum science. Screw dislocations, or vortices, are 
a common dislocation type. Wave packets with nested vortices find applications in fields as 
diverse as cosmology, biosciences, or solid state physics B [3L f|, As striking examples, 
they are at the heart of schemes to generate engineered quNits in quantum information sys- 
tems in higher dimensional Hilbert spaces || 0] , are believed to be essential for the onset of 
superfluidity in Bose-Einstein condensates (BEC) [|[ ||, [H] , or allow tracking the motion 



of a single atom [[12 



A recent study of the motion of vortex lines governed by both linear and nonlinear 
Scrodinger equations describing the dynamics of atoms in harmonic traps revealed that the 
topological features of vortex dynamics are to large extent universal ||13|| . The dynamics of 
the vortices nested on localized wave packets depends on the evolution of the host beam, 
and on the interferences and interactions between the vortices ]14| . Multiple vortices nested 
on the same host typically follow dynamical evolutions which might include large vortex- 
drifts that destroy their initial arrangement, and vortex-pair annihilations that destroy the 
vortices themselves. Vortex evolutions are particularly complex in strongly nonlinear media, 
such as BEC, where the vortices can interact with each other. Therefore, a fundamental 
question arises about whether stationary or quasi-stationary vortex clusters or lattices JTSJ . 
made of vortices with equal and with opposite topological charges exist. To isolate the pure 
vortex features from the dynamics solely induced by the evolution of the host wave packet, 
it is convenient to study wave fields confined by suitable potentials, as in weakly-interacting 
trapped BEC. 

In this paper we show that vortex clusters with multipolar global wave fronts nested in 
wave fields confined by trapping potentials are non-stationary, when the number of vortices 
and their location are not constant during dynamical evolution, or at best flipping, when the 
vortices periodically flip their topological charges through extremely sharp Berry trajectories 
T6fl . In the former case, multiple vortex revivals mediated by Freund stationary point 



bundles ||17|| , that carry the necessary Poincare-Hopf indices can occur. In contrast, 
we find that a rich variety of fully stationary vortex clusters made of an increasing number 
of vortices do exist. The important point we put forward is that these clusters are globally 
linked, rather than products of independent vortices. Also, they feature a monopolar global 
wave front. We also show that the clusters exist and are robust in nonlinear systems such 
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as interacting BEC. 

We thus address the slowly-varying evolution of generic wavef unctions governed by the 
paraxial wave equation 

iA z = LA + Af(A), (1) 

where A is a complex field, £ is a two-dimensional linear differential operator containing 
a trapping potential and M{A) takes care of any nonlinear contribution. We assume the 
trapping potential to be harmonic thus C = —1/2 {d 2 . + d 2 ) + {n x x 2 + n y y 2 ) . To be specific, 
when M{A) ~ |^4| 2 A, this equation models the propagation of a light beam guided in a Kerr 
nonlinear graded-index medium and the mean-field evolution of a two-dimensional trapped 
BEC at zero-temperature (where n xy are proportional to the trap frequencies in appropriate 
units). Here we will consider only the symmetric case; hence n x = n y = 2. For convenience, 
from now on we will split A(x, y; z) = F(x, y; z)V(x, y; z), taking the host packet F as given 
by the fundamental mode of the trapping potential F(x,y;z) = exp(— x 2 — y 2 ) exp(— 2iz). 
In the linear case, the function V(x,y;z) carries all the essential information about the 
solutions and in particular about vortex dynamics. Here we will consider the evolution of 
polynomial initial data for V corresponding to (multi) vortex solutions of Eq. (fj). Such 
solutions can be expressed as finite series and, as will be clear later, all of them must be 
periodic or stationary. 

Let us first consider the linear (Af(A) = 0) evolution of vortex-clusters built as products 
of n independent single-charge vortices: V(x, y; z = 0) = nfc=i[ x ~ x k + ick(y — Vk)}, where 
{xki Vk) are the locations in the vortex cores in the transverse plane, and = ±1. None of 
the above product-vortex clusters is found to be dynamically stationary. On the contrary, 
the number of vortices and their location is found to vary during evolution, so that the initial 
vortex structure is destroyed. These results can be illustrated by examining the evolution 
of the 4- vortex cluster: V(z = 0) = (x + a + iy) (x — a + iy)(x — iy — ia)(x — iy + ia), which 
contains two vortices with positive topological charge and two vortices with negative charge 
in a symmetrical geometry. This cluster features a quadrupolar global wave front, as is 
revealed by calculating the gradient of the wave front $ far from the cluster core, to obtain 
|V$| ~ 1/p 3 , where p is the polar coordinate, similar to the corresponding electrostatic 
multipole fl9fl . Substitution into flU) yields 

V(x,y;z) = [(x 2 + y 2 )(x 2 + y 2 + 2e 4iz -2)}e- 8iz 
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+Aia 2 xye~ Aiz + 1/2(1 - e^ iz f - a 4 . (2) 

One finds three different regimes of evolution, as shown in Fig. 1: Vortex drifts, and vortex- 
pair annihilations and revivals take place, so that depending on the value of the geometrical 
parameter a, the total number of vortices, n, hosted in the wave field during propagation 
can oscillate between: (i) 4 and 8 [see Fig. 1(a)]; (ii) 4, and 8 [see Fig. 1(b)]; (iii) 4 and 
[see Fig. 1(c)]. When n — 0, it is understood that all vortices have annihilated each 
other. Analogous evolutions were found for octupoles and higher-order multipoles. Only 
dipoles can be made quasi-stationary, but flipping, when the corresponding vortex twins 
periodically flip their topological charges. Thus, the main conclusion reached is that the 
interference between the constituent vortices of all the product-clusters produces beatings 
between the normal modes of the system, rendering the clusters non-stationary. 

The key insight we put forward in this paper is that such beatings are not associated to 
the intrinsic or local properties of the individual vortices, but to the very way the vortices 
are globally linked in the host wave packet. As an example, let us consider the evolution of 

V(x, y; z = 0) = x 2 + y 2 - a 2 + 2ixy, (3) 

which contains four vortices located at the same positions and having the same charges 
as those of the vortex-quadrupole considered above [see Fig. 2]. However, in this cluster 
the vortices are intimately linked to each other, rather than individually nested in the host 
F. This fact manifests itself in the global wave front of the cluster, which behaves as 
|V$| ~ cos(20)/p + 0(1/ p 3 ), and thus features a monopolar decay almost everywhere. In 
this case, the vortex evolution is given by 

V(x, y; z) = (x 2 + y 2 - 1/2 + 2ixy)e~ Aiz + 1/2 - a 2 , (4) 

an evolution that when the cluster is constructed with a = 1/ \/2 does become fully station- 
ary. 

The above 4- vortex cluster is not an isolated case, but rather an example of whole existing 
families of fully stationary vortex structures made of globally linked vortices. In fact, the 
solutions of Eq. ([I]) with AT (A) = have the form: 

oo 

A(x, y;z)=J2 CMH^xV^my^e-^-y'e- 2 ^ 1 ^, (5) 

k,l=0 
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where Hj are the Hermite polynomials. Therefore, the evolution of initial data of the form 
V(x, y; z = 0) = Ylk=o C k H k (^)H n . k (ri) for any C k G C and £ = x^/2, 77 = yV2, is given by 



On physical grounds, this simple mathematical result shows that all the stationary clusters 
are made of globally-linked vortices. Equation @ allows us to build a variety of structures, 
to be termed Hermite, or H-clusters, whose key features we discuss in what follows. 

Perhaps the simplest type of H-clusters are those with a n x n matrix geometry thus 
containing n 2 vortices. These clusters can be generated by using as initial data, for example, 
V nxn (x, y; z = 0) = H n (£) + iH n (r}). In this function the vortex charges alternate throughout 
the matrix and the vortex locations are dictated by the zeroes of the particular Hermite 
polynomials involved. In general these vortex matrices are not regular, the distance between 
vortices varying along the matrix. However, in the particular cases with n = 2 and n = 3 
the matrix is regular. The n = 3 case is shown in Fig. 2(a). Notice that a 2 x 2 matrix 
cluster can be generated either with V(x,y;z = 0) = H 2 (£) + iH 2 (r]), or with V(x, y; z = 
0) = H 2 (£,) + H2{rj) + iHi(£)Hi{rj). Actually, this latter possibility generates the stationary 
4- vortex cluster discussed earlier [see Fig. 2(b)]. 

One can also build m x n (m ^ n) stationary vortex matrices. A possible choice 
for the vortex function generating such a vortex-matrix is V mxn (x, y; z = 0) = H m (£) + 
iH\ m _ n \(^)H n (j]). As an example we show the 4x2 vortex matrix in Fig. 2(c). In contrast 
to the n x n matrices, in the general case the topological charges carried by the vortices of 
the m x n matrices do not alternate sign throughout the matrix. An important subclass of 
the m x n vortex matrices are the m x 1 cases, to be termed vortex arrays. They consist 
in m co-linearly displaced vortices of the same topological charge. Figures 2(d) and 2(e) 
show illustrative examples. The simplest array is the vortex- twin shown in Fig. 2(d): A 
pair of identical vortices that, contrary to the vortex dipole which either undergoes periodic 
annihilations and revivals or charge flip-flops, can be made fully stationary. 

The n x n matrices are either chargeless for even n, or carry a single net charge for odd 
values of n, while the mxl arrays carry a m total topological charge. In any case, the wave 
front of all the H-clusters is found to feature a monopolar decay (~ almost everywhere. 

More complex H-clusters also exist, and a full classification of all the possibilities falls 
beyond the scope of this paper. However, an example of one of such exotic H-clusters 



n 




(6) 
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is displayed in Fig. 2(f), which corresponds to the cluster built with V(x, y;z = 0) = 
H 3 (r/) + i(if 3 (£) + H\ (77)^2(0) • The rich variety of possibilities contained in Eq. (||) is 
clearly apparent. 

An interesting issue is the existence and stability of vortex-clusters in the presence of 
nonlinear cubic interactions such as those appearing in the propagation of beams in Kerr 
media or in the dynamics of BEC. To ease the comparison with BEC literature we choose now 
n x = n y = 1/2 and N(A) = U\A\ 2 A p0| . In this context the evolution variable is denoted by 
t instead of z. With this choice of parameters the range of U values experimentally accesible 
for the two-dimensional case is < U < 10 2 — 10 3 |2l] . 

We have studied several particular examples to verify that these structures indeed exist 
and are stable in the nonlinear regime. We have taken as initial data several linear configu- 
rations such as a single vortex, dipole systems and the 4- vortex cluster given by Eq. (^) and 
evolved them for UN = 10 (N = J \A\ 2 dit is the wave function norm) using an standard 
split-step integrator. It is found that, although the background performs oscillations and 
the vortex locations oscillate around their equilibrium positions, the vortex clusters remain 
stable (see Fig. 3). We have also searched for stationary solutions of Eq. ([!]), of the form 
A(x,y;t) = e lXt ip(x,y). To do so we have used a steepest descent method to minimize the 
functional P^fl : 

pm JV(-A-A + r 2 + f/M 2 )^x 
m j>| 2 dx ' {) 

whose minima (except for ip = 0) coincide with the stationary solutions of Eq. ([!]) for a 

given value of A. For instance, taking as initial data for the minimization process the linear 

4-vortex cluster and setting A = 8.0 and U = 100, we found a stationary 4-vortex cluster 

solution (see Fig. 4) with norm N — J |?/>| 2 <ix ~ 1.6005 (thus the product UN ~ 160 which 

lies into the fully nonlinear regime). We have verified that this solution is robust under time 

evolution when small perturbations are added. These evidences show that the existence of 

H-clusters in a BEC should be experimentally accessible, at least from the dynamical point 

of view. 

To conclude, we stress that the constituent vortices of the H-cluster are globally linked, 
rather than products of independent vortices. Following this idea it is possible to generate 
a variety of additional novel structures with fascinating properties. A nice new example are 
the circular vortex necklaces generated at the intersection between the circle x 2 + y 2 — a 2 = 0, 
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where dte(V) = 0, and the lines y ± tan(2k7c/n)x = 0, k G N, where Qm(V) = 0. Those 
are quasi-stationary, purely flipping clusters made of n vortices (see Fig. 5), a feature so far 
only known to occur with vortex-dipoles. Once again, the vortices forming the necklace are 
intimately linked and do not exhibit a n-polar wave front. The exploitation of such intrinsic 
linking might open new opportunities in classical and quantum systems based on topological 
light and matter waves. The first challenge is the demonstration of the generation of the 
clusters, by suitable computer-generated holograms p3j in Optics and phase- imprinting 
techniques in BEC [|1 . 
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Figure captions 

Fig. 1: Evolution of a vortex quadrupole, constructed as the product of four single-charge 
vortices. Upper row: Number of vortices n as a function of z, for different values of the 
initial quadrupole size, (a) a = 0.5, (b) a = 1.1 and (c) a = 1.2. Bottom row: Intensity 
snapshots corresponding to points labeled A, B and C in (c). Black-filled circles: positive 
vortices; white-filled circles: negative vortices. 

Fig. 2: H-cluster zoology. Shown are several examples of stationary vortex matrices and 
vortex arrays that can be constructed. See text for details. Lines show zero crossings of 
3?e(V) (full lines) and Qm(V) (dashed). Features as in Fig. 1. 

Fig. 3: Stable evolution of initial data given by Eq. ([|) for UN = 10. Upper row: Intensity 
plots; bottom row: Interference fringes. Spatial region spanned is [—4,4] x [—4,4]. 

Fig. 4: Linear 4-vortex cluster [Eq. @] vs. its nonlinear stationary version for U = 
100, A = 8. (a) Plots of \ip(x, y = 0)| 2 for the linear (dashed line) and nonlinear (solid line) 
cases. (b,c) Surface plots of 1^(^,2/) | 2 for (b) the linear and (c) nonlinear situations. The 
vortex locations and topological charges are indicated by plus and minus signs. 

Fig. 5: Evolution of a flipping n = 8 circular vortex- necklace. Vortex pattern at (a) z = 
and (b) z = 37r/16. 
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